Interference and zero-bias anomaly in tunneling between Luttinger-liquid wires 
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We present theoretical calculations and experimental measurements which reveal the Luttinger- 
liquid (LL) nature of elementary excitations in a system consisting of two quantum wires connected 
by a long narrow tunnel junction at the edge of a GaAs/AlGaAs bilayer heterostructure. The 
boundaries of the wires are important and lead to a characteristic interference pattern in measure- 
ments on short junctions. We show that the experimentally observed modulation of the conductance 
oscillation amplitude as a function of the voltage bias can be accounted for by spin-charge separa- 
tion of the elementary excitations in the interacting wires. Furthermore, boundaries affect the LL 
exponents of the voltage and temperature dependence of the tunneling conductance at low energies. 
We show that the measured temperature dependence of the conductance zero-bias dip as well as 
the voltage modulation of the conductance oscillation pattern can be used to extract the electron 
interaction parameters in the wires. 

PACS numbers: 73.21.Hb,71.10.Pm,73.23.Ad,73.50. Jt 
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I. INTRODUCTION 

Quasi-one-dimensional (ID) structures with gapless 
electronic excitations, such as carbon nanotubes, quan- 
tum HaU edge states, and confined states at the edge of a 
quantum well heterostructure (i.e., quantum wires), pos- 
sess unique properties which cannot be described by Lan- 
dau's Fermi-liquid theory. Even small electron-electron 
interactions in a ID confinement make inadequate the 
picture based on the existence of long-lived fermionic 
quasiparticles which can be mapped onto single-particle 
states in a free-electron gas. A powerful framework for 
understanding universal properties of ID electron sys- 
tems was put forward by the formulation of Luttinger- 
liquid (LL) theoryi (For a review see Refili) The spec- 
tral density, A{k,Lu), of the one-electron Green function 
in a Luttinger liquid is fundamentally different from that 
of a Fermi liquid: While the latter has one quasiparticle 
peak, the former has two singular peaks corresponding 
to the charge- and spin-density excitation modest^ 

Tunnel-coupled quantum wires of high quality created 
at a cleaved edge of GaAs/AlGaAs double-quantum- well 
heterostructures appear to be an exceptional tool for 
probing spectral characteristics of a ID systemiSiSii It 
is achieved^ by measuring the differential conductance 
G{V, B) as a function of the voltage bias between the 
wires, V, and magnetic field oriented perpendicular to 
the plane of the cleaved edge, B, allowing for simultane- 
ous control of the energy and momentum of the tunnel- 
ing electrons. In a recent article- we demonstrated that 
the picture of noninteracting electrons can be used with 
great success to explain some of the most pronounced fea- 
tures of the conductance interference pattern arising from 
the finite size of the tunneling region. Taking electron- 
electron interactions into account was shown to explain 
experimentally observed long-period oscillation modula- 



tions in the V direction, which can be understood as 
a moire pattern arising from spin-charge separation of 
electronic excitations. In this paper we use LL formal- 
ism to further investigate an interplay between electron 
correlations and the finite length of the tunnel junction, 
which allows us to understand peculiarities of the oscil- 
lations and the zero-bias anomaly in the measured tun- 
nehng conductance G{V,B). 



II. EXPERIMENTAL METHOD 

In this section we describe the means by which we mea- 
sure the tunneling conductance through a single isolated 
junction between two parallel wires. 



A. Fabrication of the samples 

The two parallel ID wires are fabricated by cleaved- 
edge overgrowth (CEO), see Fig.^and Ref.0. Initially, a 
GaAs/AlGaAs heterostructure with two closely situated 
parallel quantum wells is grown. The upper quantum well 
is 20 nm wide, the lower one is 30 nm wide and they are 
separated by a 6 nm AlGaAs barrier about 300 meV high. 
We use a modulation doping sequence that renders only 
the upper quantum well occupied by a two-dimensional 
electron gas (2DEG) with a density n w 2 x 10^"'^ cm~^ 
and mobility w 3 x 10^ cm^V^^s^^. After cleaving 
the sample in the molecular beam epitaxy growth cham- 
ber and growing a second modulation doping sequence, 
two parallel quantum wires are formed in the quantum 
wells along the whole edge of the sample. Both wires 
are tightly confined on three sides by atomically smooth 
planes and on the fourth side by the triangular potential 
formed at the cleaved edge. 
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FIG. 1: Illustration of the sample and the contacting scheme. 
The sample is fabricated using the CEO method. The parallel 
ID wires span along the whole cleaved edge (right facet in the 
schematic). The upper wire (UW) overlaps the 2DEG, while 
the lower wire (LW) is separated from them by a thin AlGaAs 
barrier (AlGaAs is dark gray in the schematic, electron liquids 
are light gray). Contacts to the wires are made through the 
2DEG. Several tungsten top gates can be biased to deplete 
the electrons under them: We show only Qi, here biased to 
deplete the 2DEG and both wires, and Q2, here biased to 
deplete only the 2DEG and the upper wire. The magnetic 
field B is perpendicular to the plane defined by the wires. The 
depicted configuration allows the study of the conductance of 
a tunnel junction between a section of length L of the upper 
wire and a semi-infinite lower wire. 



Spanning across the sample are several tungsten top 
gates of width 2 /zm that lie 2 /zm from each other (two 
of these are depicted in Fig. P). The differential conduc- 
tance G of the wires is measured through indium contacts 
to the 2DEG straddling tungsten top gates. While mon- 
itoring G with standard lock-in techniques (we use an 
excitation of 10 /uV at 14 Hz) at T = 0.25 K, we decrease 
the density of the electrons under the gate by decreasing 
the voltage on it (Vg). At Vg = V2B, the 2DEG depletes 
and G drops sharply, because the electrons have to scat- 
ter into the wires in order to pass under the gate. For 
V2D > Vg > Vu the conductance drops stepwise each time 
a mode in the upper wire is depleted.^*' In this voltage 
range, the contribution of the lower wire to G is negligi- 
ble because it is separated from the upper quantum well 
by a tunnel barrier. When Vg = Vjj, the upper wire de- 
pletes and only the lower wire can carry electrons under 
the gate. This last conduction channel finally depletes at 
Vl and G is suppressed to zero. 



B. Measurement on an isolated tunnel junction 

The measurements are performed in the configuration 
depicted in Fig. ^ The source is the 2DEG between 



two gates, Qi and Q2 in Fig. the voltages on which 
are Vi < Vl and Vl <V2 <Vu, respectively. The upper 
wire between these gates is at electrochemical equilibrium 
with the source 2DEG. This side of the circuit is sepa- 
rated by the tunnel junction we wish to study from the 
drain. The drain is the 2DEG to the right of Q2 (the semi- 
infinite 2DEG in Fig. ^ and it is in equilibrium with the 
right, semi-infinite, upper wire and with the whole semi- 
infinite lower wire in Fig.^ Thus, any voltage difference 
{V) induced between the source and the drain drops on 
the narrow tunnel junction between the gates. This con- 
figuration gives us control over both the energy and the 
momentum of the tunneling electrons, as explained be- 
low. An additional gate lying between Qi and Q2 (not 
shown in Fig. ^ allows us to deplete the 2DEG in the 
center of the source, thus reducing the screening of the 
interactions in the wires by the 2DEG. 

The energy of the electrons tunneling between the 
wires is given by eV ^ — e being the electron charge. The 
tunneling process occurs along the whole length L of the 
tunnel junction. Therefore, momentum is conserved to 
within an uncertainty of order 2tt/L <C fcp, where fcp is 
a typical Fermi wave vector in the wires. We can shift 
the momentum of the tunneling electrons with a mag- 
netic field {B) perpendicular to the plane defined by the 
wires. The value of the shift is fiqB = ei?d, where d is 
the center-to-center distance between the wires. 



III. DESCRIPTION OF THE EXPERIMENTAL 
RESULTS 

In the experiment we measure the nonlinear differ- 
ential tunneling conductance G{V^ B) through a junc- 
tion between two parallel wires. The sample we re- 
port here contains four top gates allowing us to vary 
the length of the junction L by choosing different com- 
binations of gates. We have studied in detail junctions 
with L = 2, 4, 6, 10 /um as well as symmetric junctions 
(L = 00). The results presented here are from junctions 
with L — 2, 6, 10 /zm. Many of the effects that we mea- 
sure rely on the smallness of 1/L, while others (which we 
address here in detail) are present only when L is finite. 



A. Dispersions of elementary excitations in the 
wires 

By mapping out G{V^ B) we determine the dispersion 
curves of the wires. ^ These are given by the curves that 
are traced by the main peaks as seen in Fig.[2| We can un- 
derstand their gross features employing a nonintcracting 
electron picture:^ The peaks result from tunneling be- 
tween a Fermi point in one wire and a mode in the other 
wire. Since each occupied mode has two Fermi points, 
two copies of the dispersion show up in the G(V, B) scan. 
All in all, for each pair of occupied modes in the two wires 
we expect to observe four dispersions, because there are 
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FIG. 2: Plot of G(V, B) for a 10 /im junction. Higher values 
of the conductance are depicted in lighter shades: The top 
bar gives the key. 



four Fermi points involved: ±fcp„ and ±fcpj. (Indices i 
and j label various modes in the wires, u and ^~the upper 
and the lower wires.) In reality, we observe only some of 
the transitions: For example, by carefully studying Fig. [21 
one can distinguish dispersions of three modes in the up- 
per wire and five in the lower one, but only the following 
transitions seem to have a sizable signal: <-> 
1^3) ^ 1^2), and 1^2) 1^3,4,5), where the order in the 
list is of decreasing B^'' (see below); \uz) is the 2DEG 
occupying the upper quantum well.^ Such selection rules 
are related to the shape of the wave functions in the di- 
rection perpendicular to the cleaved edge. In identical 
wires, one would expect only the transitions |m„) |Z„) 
to show up, due to the orthogonality between different 
modes. In dissimilar wires, the selection rules are differ- 
ent and less strict so other transitions are observed. 

The dispersions allow us to extract the densities of elec- 
trons in each mode, '^■^(j) = {'^I'^)^fu{i)^ follows. Tun- 
neling amongst each pair of occupied modes is enhanced 



near 1^ = at two values of i? > 0, where the two curves 
in G(V, B) cross. In the first (in the following referred to 
as the "lower crossing point"), which occurs at BW the 
direction in which the electrons propagate is conserved 
in the tunneling process. In the second (referred to as 
the "upper crossing point" ) , the Lorentz force exerted by 
B2'' exactly compensates for the momentum mismatch 
between oppositely moving electrons and the direction of 
propagation of the tunneling electrons reverses. In wires 
with vanishing cross section, these crossing points occur 
at 
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In principle Eq. can be used to extract the densities 
of the modes, regardless of electron-electron interactions 
in the wires^ or mesoscopic chargingii^ that can merely 
smear them at a finite voltage bias. In realistic wires that 
have a finite cross section, finding the densities is ham- 
pered by the weak magnetic field dependence that they 
acquire. This difficulty is overcome by a simple fitting 
procedure that we have developed: We assume that all 
the modes in a wire have the same field dependence, a 
reasonable assumption for our tight-confining potential 
in the growth direction of the quantum wells. We then 
guess the B = occupations of the modes in each wire, 
n*^j(0) and (0), and calculate their field dependences. If 
the resulting dispersions do not cross at B^^^ , we adjust 

n^(0) and nf{0) and repeat the procedure. This is done 
iteratively for all the crossing points that we see, because 
changing the occupation of one mode affects the field de- 
pendence of all the other occupations in a wire. The 
dispersion that we use is that of noninteracting electrons 
in a finite quantum well, in the presence of an in-plane 
magnetic field. Such a dispersion depends only on the 
width and depth of the well and on the band mass of 
electrons in GaAs. 

In every case we have studied, we see clear deviations 
of the measured dispersions from the calculated nonin- 
teracting ones at a finite bias. In particular, we find that 
the velocities of some excitations are enhanced relative 
to the Fermi velocities Wfu(/)- The former are given by 



1 dV 



(2) 



(along the observed main peaks), while the latter can 
be obtained by the calculated slope of the (noninteract- 
ing) dispersions at the Fermi points. This velocity en- 
hancement is thought to correspond to the charge-density 
modes and can be accounted for by electron-electron in- 
teractions in the wiresi^ 

The ability to determine the dispersion relations re- 
lies on the high quality of the junctions to sustain 
momentum-conserving tunneling. Momentum relaxation 
ensues as soon as invariance to translations is broken. 
The most obvious mechanism by which this happens is 
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FIG. 3: Nonlinear conductance oscillations at low field from 
a 2 nm junction, (a) shows the oscillations as a function of 
both B and V. (A smoothed background has been subtracted 
to emphasize the oscillations.) The brightest (and darkest) 
lines, corresponding to tunneling between the lowest modes, 
break the V-B plain into regions I, II, and III. Additional 
positively-sloped bright and dark lines in II arise from other 
ID channels in the wires and axe disregarded in our theoretical 
analysis. Also present is a slow modulation of the strength 
of the oscillations along the abscissa, (b) Absolute value of 
the peak of the Fourier transform of S'^'^^'^G (V, 5^+^ at 
a fixed V in region II as a function of V. (See Sec. lIVB II for 
definition of S, /3 and other details.) Its slow modulation as 
a function of V is easily discerned. 



the finiteness of L. We find that we indeed observe its 
effects. The second mechanism is the disorder inherent 
to all semiconductor devices, some effects of which seem 
to also be observed. 



B. Oscillations 

The most spectacular manifestation of the breaking of 
translational invariance is the appearance of a regular 
pattern of oscillations away from the dispersion curves. 
Figs.|21i and^Ji are typical examples of the patterns that 
we measure at low magnetic field. In this range of field, 
the lines that correspond to the dispersion curves appear 
as the pronounced peaks that extend diagonally across 
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FIG. 4: Same as Fig. |21but for a 6 nm junction. Note that 
the oscillations are approximately three times faster than in 
Fig. 13 as expected from Eq. . For this junction, there are 
several additional side lobes present on the left of the principal 
peaks, unlike in the case of the shorter junction in Fig. |3 



the figures. In addition to these we observe numerous 
secondary peaks running parallel to the main dispersion 
curves. These side lobes always appear to the right of 
the wire dispersions, in the region that corresponds to 
momentum conserving tunneling for an upper wire with 
a reduced density. As a result, we see a checkerboard 
pattern of oscillations in region I, a hatched pattern in 
region II, and no regular pattern in region III (see Figs.|2K 
and^ for the definitions). 

The interference pattern also appears near the upper 
crossing point at high magnetic field. A typical example 
is shown in Fig. \^ 

The frequency of the oscillations depends on L. When 
L is increased from 2 /xm. Fig. [SJ to 6 /xm, Fig. ^ the 
frequency in bias (AF) and in field {AB) increases by 
about a factor of three. The period is approximately 
related to the length of the junction through the formula 



AVL/vF ^ ABLd = 



(3) 



where = ^Trh/e is the quantum of flux. 

A close examination of the low-field oscillations reveals 
an interesting behavior of their envelope. Notable is the 
suppression of G(V, B) near V = which is independent 
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FIG. 5: G{V, B) near the upper crossing point for a 6 fim 
junction. In this measurement, a central 2 jj,m gate midway 
between C/i and Q2 is biased to deplete all upper wire modes 
except the lowest one. One can see a pattern of oscillations 
around the dispersion peaks. 

of field. Also visible are faint vertical gray stripes, where 
the amplitude of the oscillations in the B direction is 
reduced. The modulation of the oscillation amplitude, 
as a function of V, is shown in panels (b) of figures|3|and 
0] The oscillatory part of G thus depends on V on two 
major scales: The faster scale (0.5 mV for L = 2 /im) 
corresponds to the oscillations described by Eq. (PJ . The 
slower scale (2 mV for L = 2 /im) governs the distance 
between the stripes of suppressed G{V, B) parallel to the 
field axis, including the zero-bias suppression. Like the 
fast scale, the slow scale is roughly inversely proportional 
to the lithographic length of the tunneling region. 



C. A dip in the tunneling conductance 

Prominent in all scans that have high enough resolu- 
tion in is a strong suppression of the conductance near 
1^ = at all magnetic fields. The width of this conduc- 
tance dip is of order of 0.1 mV, see Figs.|21and|Sl The size 
of the dip is very sensitive to temperature, as depicted in 
Fig. El and it exists for T < 1.0 K. 



IV. THEORY AND DISCUSSION 

The ID modes in the upper quantum well are cou- 
pled to the 2DEG via an elastic 1D-2D scattering which 
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FIG. 6: Zero- voltage dip of the tunneling conductance G as a 
function of temperature on a log-log scale. The circles show 
measurements on a 6 fim junction at B = 2.5 T, the lines 
are a fit using G oc T" for V — 0. The dashed line is the 
result for a = abuik(5i) = 0.07 while the solid line is the re- 
sult for a — Qcnd(5i) ~ 0.35, with gi = 0.59 (and gu = 1) 
in Eqs. I|49|l and 1)50^ . respectively; see Sec. liv e 21 for a dis- 
cussion. Insets: G{V) for T = 0.24 K and T = 0.54 K (the 
temperature dependence was generated from the V — point 
of such scans). The curves were calculated with Eq. I51II 
and using the above value of gi extracted from the fit of 
the temperature dependence of the dip. [We obtained Fa{x) 
by convoluting the derivative of the Fermi distribution in 
the 2D leads, [l/(4fcsr)]sech^[el//(2fcsr)], with the finite- 
temperature tunneling density of states in the lower wire, see 
Eq. (5) in Ref. [T3 .] The dashed lines correspond to the ctbuik 
value of the exponent while the solid lines to acnd- 



ensures a good electronic transfer between the extended 
and confined states of the wellii^ In addition to tunneling 
between the confined states in the wires, if the extended 
states have an appreciable weight at the edge, there will 
be a direct transition from the 2DEG to the lower wire. 
With this in mind, wc separate the total current into two 
contributions, one due to tunneling between ID bands 
and the other due to direct tunneling from the 2DEG. 
As explained in Sec. IIII Al each of the wires carries sev- 
eral ID modes. In our analysis and comparison with the 
experiment, we will only consider the transition between 
the lowest ID bands of the wires (i.e., the bands with the 



6 



U,(x) 



U (x) 




FIG. 7: Schematic picture of the theoretical model. The 
upper wire is formed by a potential well Uu{x) created by 
gates Qi and Q2 (shown in Fig. and the lower wire is semi- 
infinite with the left boundary Ui{x) at gate Qi. 'ip{x) is an 
electron wave function in the upper wire. The energy and 
momentum of the tunneling electrons are governed by the 
voltage bias V and magnetic field B. 



\li), and the direct 
2), which both have 



largest Fermi momentum), 
tunneling from the 2DEG, jua) 
a strong signal, as seen in Fig. |21 In each wire, the ID 
modes interact with each other, but since the bands have 
very different Fermi velocities, we treat them indepen- 
dently. This is a reasonable approximation, as explained 
in Appendix 1X1 

The geometry for our theoretical description is shown 
in Fig. [71 The potentials Uu{x) and Ui{x) are felt by 
electrons in the upper and lower quantum wires, respec- 
tively. The electrons in the upper wire are confined to a 
region of finite length by potential gates at both its ends 
(see the source region in Fig.^l. One of these gates (Gi) 
causes the electrons in the lower wire to be reflected at 
one end, but the other {Q2) allows them to pass freely 
under it. The effective tunneling region is determined by 
the length of the upper wire, which is approximately the 
region |x| < L/2 in Fig. [3 The magnetic field, B, gives 
a momentum boost hqs — eBd along the i-axis for the 
electrons tunneling from the upper to the lower wire. 

First, we develop a general formalism in Sec. IIV Al 
We then apply it to study the conductance interference 
pattern in Sec. IIV 61 and the zero-bias anomaly regime in 
Sec. Hvn 



A. General formalism 

Let us first consider transport between two ID bands 
in the wires. We use the following model Hamiltonian to 
study the intermode tunneling in the system: 



H = 



V — U.l 



E 



H 



1D-2D 



(4) 



Hq is the kinetic energy of the electrons [v = u (l) la- 
bels the upper (lower) wire], H^^^ [-^hit] describes spin- 
independent electron-electron interactions in (between) 
the wires, i?iD-2D is an effective Hamiltonian for the ID- 
2D scattering of electrons in the top quantum well, and 
^^tun is the tunneling Hamiltonian: 



^int 



WFi/E / 

^E j dkV^u'{k) [2pRsu{k)pLs'u'[~k) 

ss' 

+ PRsi>{k)pRs'u'[~'k) + PLsu{k)pLs'u'{-k)] , 

AE/ ^2;*L*sie 



(5) 



(6) 
(7) 



where s and s' are spin indices, sv is the spin-s elec- 
tron field operator, rsv and lsv are the field oper- 
ators for the right and left movers, respectively, sv = 

is the density-fluctuation operator for the spin-s right 
movers (and analogously for the left movers), and 
Vuv'{k) — J dxe^^^V^u'ix) is the Fourier transform of the 
two-particle interaction potential V^yiix). Writing 7?int 
in terms of the interactions between electrons of flxed 
chirality in Eq. JHl is possible after restricting electron 
correlations to small momentum transfer scattering, e.g., 
if V{k) (X exp(— rdfel) with I/tc kp. (By making 
this approximation we disregard backward and Umklapp 
scattering processes, which are thought to be unimpor- 
tant in our cleaved-wire structure, see, e.g., Ref. (|.) 

The 1D-2D scattering randomizes the direction of the 
ID electrons in the top quantum well with a mean free 
path Z1D-2D ~ 6 fim.-^^ In infinite wires, this weak scat- 
tering can be taken into account by rounding the ID 
electron-gas spectral function by a Lorentzian of half 
width r — l/(2riD-2D), where riD_2D is the 1D-2D scat- 
tering time. 

If there were no interactions between the wires, i.e., 
Vid = 0, low-energy spin and charge excitations in each 
wire would propagate with velocities = v-pi, and Vcv — 
VFu/gv, respectively. The parameters can be obtained 
by bosonization as 



9u 



2K.(0) 



-1/2 



< 1, 



(8) 



in the case of repulsive interactions, V^,y{0) > 0. In the 
limit of a free-electron gas, V^^{0) = 0, = 1. 

We treat tunneling between the wires to lowest order in 
perturbation theory. Mesoscopic charging effects, such as 
discussed in, e.g., Ref. I 111 are disregarded in our analysis. 
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The current (for electrons of each spin) is given by 



1. Asymmetry due to soft boundaries 



I = e\X\' 



dxdx' 



(9) 



where C{x,x';t) is a two-point Green function 



x'; t) = ( U\-^u{x, t), *t 0) 



(10) 



In the hmit of vanishing interactions between the wires, 
it reduces to 

C{x,x';t) = G>{x,t;x',0)Gf{x',0]x,t) 

-G<ix,t;x',0)G>{x',0;x,t) (11) 

expressed in terms of the one-particle correlation func- 
tions 

G>{x,t;x',t') = -t{'f,{x,t)-i>l{x',t')) , (12) 
G<{x,t;x',t') = t{^lix',t')^',{x,t)) . (13) 

Note: Throughout this paper, as in the above equations, 
the correlation functions are defined for electrons with 
a fixed spin orientation and the spin index is therefore 
omitted. 

The results of this section also hold for direct 2DEG- 
ID tunneling, if we define 'i'su{x, t) as the field operator 
for the 2DEG at the edge of the upper quantum well. 



B. Interference pattern 



As discussed in Scc. lIIIBl the breaking of translational 
invariance due to the finite size of the tunneling junction 
can result in an oscillatory dependence of the conduc- 
tance G on voltage bias V and magnetic field B. In this 
section we discuss in detail this behavior that arises due 
to interference of electrons tunneling through a finite- 
sized window. We show that our theoretical framework 
can quantitatively explain the conductance oscillations 
observed near the crossing points. 

In the following, we mainly focus on the analysis of very 
distinct interference patterns measured at low magnetic 
fields (as in Figs-ElandlH). In Sec UVBH we briefiv com- 
ment on the conductance near the upper crossing point 
at high fields (as in Fig.O. It appears likely that while 
in the former regime the translational invariance is bro- 
ken due to the finiteness of the tunneling region only, in 
the latter case some other mechanisms may also play a 
prominent role. 

In the actual experiments, several ID electron modes 
are occupied in the wires. Here we consider only tunnel- 
ing between modes which have the lowest energy of trans- 
verse motion, and hence the largest Fermi momentum 
along the wire, namely and These modes have 
densities that differ by only a few percent (see Ref.|3). We 
thus make a simplifying approximation vfu = ufi = '^f, 
which is justified by the measured dispersion slopes^ 



In Ref. la we showed that the observed asymmetry of 
the secondary oscillation peaks on the two sides of the 
main dispersion curves (see Figs. 01 can be explained 
within a noninteracting electron picture and assuming a 
soft confining potential Uu{x) for the upper wire. Here we 
will employ the model developed there to quantitatively 
study the form of [/„ (x) . 

Using the phenomenological tunneling Hamiltonian 
(0, we express the current through the junction at zero 
temperature 

I^sgn{V)J2\M{n,qB,V)\^ (14) 

m 

in terms of the tunneling matrix element 

M{n,qB,V)= y"dxV':(x)e^*««>fe,(x) (15) 

between the upper-wire state ipn and the lower-wire state 
ipki , the energy of which is lower by eV. The summa- 
tion in Eq. H14I) is over the integers [sgn(l/) — l]/2 < 
m < e\V\L/{nfiVF) denoting the offset of the tpn index 
n — UF + sgn{V)m with respect to the state ipn^ just 
below the Fermi energy of the upper wire (and lineariz- 
ing the dispersions near the Fermi points, assuming e\V\ 
is not too large). The current 1)14(1 can be expressed 
by a single sum because the states of the confined up- 
per wire are discrete, while the states in the lower wire 



can be indexed by a continuous wave vec- 



tor ki . (As in Ref. |3, it is assumed that the left boundary 
Ui{x) of the lower wire lies outside the tunneling region.) 
Since the Zeeman energy in GaAs is small, we ignore the 
spin degrees of freedom. 

We argued^ that for practical purposes of understand- 
ing our measurements, the sum in Eq. (|14|l can be re- 
placed by an integral 



eV 



/(X / dt\M{EFu + e,qB,V)Y 



(16) 



labeling states in the upper wire by energy e with re- 
spect to the Fermi energy E-pu- For the conductance 
obtained by differentiating the current, this approxima- 
tion will smear out the (5-functions appearing when the 
chemical potential of the upper wire crosses each discrete 
energy level. Physically, such smearing can be caused 
by 1D-2D scattering and finite temperature. But even 
at low temperatures and vanishing scattering, the re- 
sult obtained by integration [Eq. IjlGI) ] will not be far 
off from that found by summation [Eq. (|14|l ] as the domi- 
nant contribution to the oscillation pattern near the lower 
crossing point comes from differentiating the summand 
in Eq. (|14|l [or correspondingly the integrand in Eq. (|16() ]. 
as explained below. 

We linearize the dispersions about the Fermi wave vec- 
tors k-pv, so that ki is given by {ki — k-pi)v-p = {e — eV)/h. 
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The wave vector inside the upper wire similarly depends 
on energy: — k-pu)v-p — e/h. The matrix element 
squared \M{Epu + can then be written as 

a sum of contributions due to tunneling between right 
movers and between left movers, 

|M(i?F„ + e, to, V)\^ = \MiK+)\^ + |M(k_)|2 , (17) 



where k± — ku — 
Afcp = k-pu - kpi- 



ki ± QB = Afcp + eV/{hvF) ± QB and 
The tunneling matrix element 



M{k) = / dxe"'''Mx)e 



(18) 



is determined by the form of the bound-state wave func- 
tion ipu{x) at the Fermi level of the upper wire. We wrote 
the right-hand side of Eq. H17|l as an incoherent sum of 
the contributions of the two chiralities. This is an approx- 
imation we make by disregarding additional interference 
arising due to the reflection of electrons in the lower wire 
under gate Gi (i.e., by the potential Ui in Fig. [71). Taking 
the latter into account does not considerably affect our 
results. 

I-^(^±)P do not depend on energy e, and the current 
(|16ll can, therefore, be written as^ 



locV [\M{k+)\^ + |M(k_)|^ 



(19) 



The differential conductance G — dl/dV corresponding 
to the current / cx V\M{V)\'^ becomes G cx \M{V)\'^ + 
Vd\AI{V)\'^ /dV. If, for example, the oscillatory compo- 
nent of |M(V)p has the form sin(const x V), the ampli- 
tude of the second term in the conductance will be 2ttN 
times larger than the amplitude of the first term after N 
periods of oscillation. The dominant contribution to the 
oscillatory component of the conductance near the lower 
crossing point is thus 



d_ 

dV 



G(xV^ [\M{k+)\^ 



(20) 



If the upper wire confining potential Uu is smooth 
enough so that the states at the Fermi energy can 
be evaluated by the WKB approximation, the form of 
M{k) [Eq. (|18() ] can be studied both numerically and 
analytically.* In the region between the classical turning 
points. 



Ipuix) 



1 



\/kuix) 



s{x) 



(21) 



where ku{x) — fcFti[l — Uu{x)/E-pu]'^/'^ and s{x) = 
dx'[kpu — ku{x')]. In the stationary-phase approx- 
imation (SPA), M{k) is evaluated near positions x^ 
{x~^ > x~) where fc„(a;^) = fcpu — and the integrand 
in Eq. 1)18(1 has a stationary phase. In the case of a sym- 
metric potential, Uu{x) = C/„(— x), the SPA gives 



M{k) cx 



: COS KX 



s{x+) - 7r/4] , (22) 
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FIG. 8: \M{K,)f obtained using wave function -i/^u for the 
100th WKB state in the potential well [Eq. ^ with /3 = 8] 
of the upper wire. The solid line is the numerical calculation, 
the dotted line is the SPA approximation [Eq. 122ll 1 and the 
dashed line shows the result for the square-well confinement, 
for comparison. 



where Q{k) is the Heaviside step- function, the prime in 
denotes the derivative. The SPA approximation H22() 
diverges for small values of k and we have to resort to a 
numerical calculation of the integral in Eq. (|18|l .* Fig. |H| 
shows the calculated |M(k)|^. 

We study the profile of confinement Uu{x) by measur- 
ing the period Ak of the | Af (k) |^ oscillations as a function 
of K. In a square well of length L, this period is given 
by 271 /L. In a soft confinement, the interference stems 
from the oscillations of the electron wave function near 
the classical turning points, so that Ak « 27r/(.T"*" — x~). 
For a potential of the form 



Uu{x) = Eyu 



2x 

T 







(23) 



(where (3 characterizes the ratio between the total length 
of the upper wire and the extent of its boundaries) 
2x+/L « {2K/k-pY/l^ for K > assuming that Afcp < 
/cp = {kpu + fcF/)/2r' = —x^ for a symmetric 

potential) and the period is therefore given by Ak k, 
{2tt / L){kY /2k)^/ ^ . Experimentally, we extracted Ak by 
measuring the distance between oscillation zeros in re- 
gion II of the interference shown in Fig. |31 In order to 
reduce the statistical uncertainty, the conductance was 
averaged along lines of constant k^, separately for posi- 
tive and negative bias. In terms of variable S = Tik+ /ed 
(which reduces to magnetic field B at zero voltage and 
vanishing A/cp), 



AS-; 



27r;i 
edL 



hk-p 
2edS 



1//3 



(24) 



This AS* is compared with the data in Fig. |51 for sev- 
eral values of (3 (we extract kp w 1.5 x 10* m^^ us- 
ing measured electron densities^). At each f3 shown in 
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FIG. 9: Period of (faster) oscillations in region II of Fig. |3 
as a function of 5 = hK+/{ed). Circles show measurements 
at positive and negative bias and the curves are fits using 
Eq. 12411 at several values of /?. The best overall fit is reached 
at /3 « 7.67 and L/Luth ~ 1.41, where Luth ~ 2 /im. 



the figure, the distance L was found by the best (least- 
square) fit of the curve (|24(l to the measurements. The 
Uthographic length for the junction was Lmh = 2 /im 
and the width o? = 31 nm. Such fitting allows us to ex- 
tract two quantities, L/Luth — 1-45 ± 0.1 and P — 8±2, 
characterizing the extent of the ID confinement and the 
sharpness of the potential- well boundaries, respectively. 
It appears that the effective length of the upper wire [de- 
fined as the distance between the classical turning points, 
see Eq. (|23|) ] L is actually about a micron longer than the 
lithographic length. This conclusion is relatively insen- 
sitive to the fitting procedure, as AS' in Eq. H24() ap- 
proaches 2nT%/{edL) for S ^ tikp/{2ed) if the exponent 
1//3 is small. The difference between L and Luth can be 
due to significant screening of the tungsten gates (which 
are positioned 0.5 /im above the junction) by the 2DEG 
in the upper quantum well, as viewed by the upper-wire 
electronic bands. 

As a consistency check for the result of the fit in Fig. |51 
we performed an analysis of the conductance oscillations 
that takes into account the dependence on f3. Accord- 
ing to the S~^^^ scaling of the oscillations' period, see 
Eq. (|23I), and the S^^'^^'^ fall-off of the their amplitude 
[which follows from Eq. H22I) , also see Ref. 8] , the function 
5i-i//3G'(V",5i+i/^) is periodic in S^+^/l^. Fourier ana- 
lyzing it at a fixed V, and setting /3 = 8, we obtained a 
main peak, the position of which depends very weakly on 
V, which corresponds to a length oi L — 2.81 ± 0.02/xm, 
in agreement with the result of the fit. In Fig. we 
plot the absolute value of the main peak, which is seen 
to decay on a scale of a few mV. We discuss this decay 
in Sec. lIVB4l For comparison, we also Fourier analyzed 
the data in Fig.^. For that we found that one has to use 
a larger value, (3 = 21.5, in order to obtain a relatively 
voltage-independent position of the peak. This value of 
/3 is reasonable because it gives approximately the same 



boundary profile for a 6 /zm (upper) wire as /3 = 8 gives 
for a 2 /im wire,i^ Again we obtained a reasonable length 
{L = 7.3 ± 0.3 /im) that varied only weakly as a function 
of V. The height of the main peak in this case is shown 
in Fig. 0Jd where it is seen to decay on a faster scale than 
for the shorter upper wire. The ratio of the scales is 
approximately the ratio of the upper-wire lengths. 



2. Modulation due to spin-charge separation 

In the following we describe how electron-electron in- 
teractions in the wires and between them affect the os- 
cillation pattern. We show our theoretical results for 
G{V,B) near the lower crossing point of the |ui) ^ \li) 
transition and compare them to measurements on 2 /im 
and 6 /xm junctions. Figs. |21 and ^ In particular, we 
find that the difference in the velocities of the charge- 
and spin-excitation modes in the double- wire system can 
account for the observed G{V, B) suppression stripes run- 
ning parallel to the i3-axis. 

As a starting point, let us consider the case when the 
interwire interactions are vanishingly small Vui <C Vu and 
the interactions in the two wires are the same, Vuu — Viu 
so that gi ^ Qu = g, as defined in Eq. ((Sj). For positive 
voltages > 0, the current is then given by 



/ = el 



xG>{x,t;x',0)Gf{x',0;x,t) . 



(25) 



At low magnetic field, the conductance has two main 
contributions, corresponding to the two edge-state chi- 
ralities. The two contributions give bright conductance 
peaks and side lobes with opposite slopes, as described 
in Sec. IIV Bl and Ref.'§. Let us discuss tunneling between 
the right movers (current due to tunneling between the 
left movers at field B equals tunneling between the right 
movers at field —B). We assume that the electron density 
in each wire varies slowly on the length scale set by the 
respective fcp (except for unimportant regions very close 
to the boundaries). The zero-temperature Green func- 
tions entering Eq. (|25|l . in this regime, can be written 



Gu,i{x,t;x',0) 



i—^ui{x, x') ^ 

1 



{z-Vct±iO+)- 



where Wc = fp/s, 7 = (5+5 ^-2)/4, z 



(26) 



and 



Tr. IS 



a short distance cutoff (i.e., l/rc is a momentum-transfer 
cutoff in the electron-electron interactions). Here, G^ is 
the Green function (|12|l for the upper wire and G; 
is the G^ Green function (|13|) for the lower wire. The 
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function $y is defined by ^i^{x,x') = ijj,y{x)ip*{x'), in 
terms of the WKB wave functions ipui^) for right-moving 
electrons at the Fermi energy in wire i/ in a confining 
potential U^{x) which must be chosen self-consistently 
to give the correct electron density. Here we assume that 
ipu{x) and Uu{x) are given by Eqs. (|21|l and (|23|l . while 

Several additional approximations are implied in using 
Eq. H26(l to calculate the tunneling current H25|l : (1) The 
weak 1D-2D scattering is neglected, (2) The voltage is 
small enough so that one can linearize the noninteracting 
electron dispersions about the Fermi-points and use LL 
theory (i.e., we disregard the curvature), (3) t <C VcL, so 
that the discreteness of the energy levels of the upper wire 
due to electron confinement within a well of length L and 
their reflection at the boundaries does not considerably 
modify the LL Green function for an infinite wire [the 
confinement, however, is manifested in the form of the 
wave function ipuix)', effects due to the discreteness are 
discussed in Sec. IIVB II in the regime of noninteracting 
electrons, and they are believed to be small]. The last 
approximation breaks down for very low voltages (and, 
correspondingly, long times) in the regime of the zero- 
bias anomaly, which is treated separately in Sec. IIV C II 

Substituting Green functions (PHf into integral (|^ . 
we obtain for the tunneling current 

/oc 
dxdx'e'^'^^~''"'>^'=~''">i:uixWu{x')h{x - x') , 
-OO 

(27) 



using the definition 
h{z) = - / dt 



cieVt/h 



{z - vpt + zO+)(z - Vet + iO+) 



Z — Vet + iVc J \ Z + Vet — ITc 



(28) 



The integrand in Eq. (|28|l has a simple analytic form: 
it has two first-order poles at t = z/{vp + iO"'") and 
t — zjivc -\- lO"*"), and two branch cuts starting with sin- 
gularities at i = (±z-|-irc)/fc. The contour of integration 
can be deformed leaving two nonvanishing contributions: 
/i(z) — hi{z) + h2{z). The first contribution, hi(z), is 
due to integration around the poles: 



hi{z) 



{vc — wf)(-z + iO+) V'"c + 2izrc 

27rie'''^^/('"™'') 
~ {vc-VF){z + iQ+) 



rl + z^[l - (vc/vf)^] + 2izrcVc/vF 



(29) 



and the second contribution, h2{z), is due to integration 
around the branch cuts. For z > 0, for example. 



h2iz) = 2isin(77r)e-'=^'''^/(''''^) 




z/Vc 



dt 



{Vct)"^ — Z^ 



^ieVt/h 



{z - Vet - irc){z - vpt - ircVp/vc) 



(30) 



In our system we expect thati^ g ~ 0.7, so that 7 ~ 
0.03 <C 1. Therefore, since Tc ~ 30 nm (the width of 
the wires) and z < L « 2 — 6 /xm, the terms of the form 
(• • •)''' in Eq. (|29|l can be safely ignored (except for the 
regime of extremely low voltages, which will be discussed 
in Sec. IIV (TT|l . Furthermore, /i2(z) ^ hi{z), so that we 
arrive at an approximation 



Hz) 



-27ri- 



{vc - VF)iz + iO+) 



(31) 



Substituting this into Eq. (|27|) . we can now evaluate the 
current. One notices that after making the approxima- 
tion H31(l . the current (|27|l becomes the same as if there 
were no electron-electron interactions but different Fermi 
velocities in the two wires, given by vf and Vc- Using 
Eqs. 12711 and H31|) . we, finally, get for the conductance 
G = dl/dV: 



G{V, B) cx 



VF 



— |M(«F)p--|Af(«;,)|2 

Vf Vc 



, (32) 



where Kp.c = 5b + A/cp -I- eV/ {Tivf^c) and M(k) is given 
by Eq. 

If the excitation velocities in the wires are nearly the 
same, vf ~ Vc = v, we can approximate the conductance 
^ by 

G (X I-77 |M(ry, vf = \M{^,f + y A \M{nf , (33) 

where M{ri,V) = M{qB + Afcp + -qV) and 77 = e/{fiv). 
This reproduces Eq. (|19ll . One can refine the form of the 
second term on the right-hand side of Eq. (|33|l (which 
can be much larger than the first term, see Sec. IIVBT|| . 
using approximation (|22|) . when the difference between 
velocities vf and Vc becomes appreciable (which is the 
case for g « 0.7): 



G cx sm[eVx+{\lvF-l/vc)m 

1/vF-l/Vc 

s{x+)] . 



UL{x+) 

X cos 2 [kx" 



(34) 



Here, n and x~^ are defined using velocity v = 2(1 /vf + 
l/wc)~^, G stands for the second contribution to the con- 
ductance in Eq. H33|l . At low bias, G ^ linearly in 
V and the term G cx |Af(K)p governs the conductance. 
This contribution is further suppressed as V"' (at zero 
temperature) in the zero-bias anomaly regime discussed 
in Sec. llVCJll 

We can generalize the preceding discussion of this 
section to include interactions between the wires, i.e., 
Vui 7^ 0. Since the quantum wires are closely spaced, the 
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interwire interactions can be sizable. Furthermore, be- 
cause the Fermi velocities in modes and are sim- 
ilar, the excitations in the coupled wires can propagate 
with velocities quite different from those in the isolated 
wires. When we take Vui into account, the dominant part 
of Green function (|10|l becomes (assuming weak interac- 
tions, in the spirit of the preceding discussion)i 



C[x,x -jt + tO^) cx - 

{Z ~ VFi,t)2(z - VFlt)2 

1 



{z - Vc-t)i+'^'- {z - Vc+t)i-'^'- 



where 



Vc± 



Vci ^ Vum 

■kTi 



(35) 



(36) 



Here, r = tt{vcu - Wc/)/2K;(0) and 0r = l/{2Vl + r^) is 
finite for nonvanishing interactions between the wires, Vcn 
are the charge-excitation velocities in the isolated wires. 
[Note that there appears to be a sign error in Ref. ^ in 
the expression for the velocities Vc± in the physical case 
of repulsive interactions Vui{0) > 0.] 

For a symmetric double-wire system, vpu = v-pi = 
V¥, Vuu = Vu, and Vcu = VcU so that r = and 
9r = 1/2. (In this case, Vc+ and Vc- become the ve- 
locities of the symmetric and antisymmetric charge ex- 
citations, respectively.) Green function then reduces to 
C oc — <i>„$;*(z — vpt)~^{z — Vc-t)~^ and we reproduce 
our main result of this section, Eq. H32|l . after replac- 
ing Vc with the antisymmetric charge-excitation velocity 
Vc~- This is natural as tunneling in a symmetric biwire 
can only excite the antisymmetric modes at low magnetic 
fields. 

In addition to the structure studied in Sec. lIVBTI for 
the system of noninteracting electrons, we now show that 
the electron-electron interactions in the wires lead to 
a modulation of the conductance oscillations along the 
voltage axis.* This modulation suppresses the contribu- 
tion G [Eq. (|34|l ] to zero in stripes parallel to the field 
axis. The distance between them is: 



AKiod = 



■nTiVc-V-p 
ex'^{vc- — v-p) 



The ratio between AVmod and the period 



Ay 



(37) 



(38) 



due to the wave-function oscillations near the turning 
points [compare to Eq. ||2J)] 



AK 



mod 



1 Vc- + Vy 1 1 + (?- 



Ay 



2 Vc 



Up 



2 1 



5- 



(39) 



can be used as an independent measure of the interaction 
parameter g_ — vf/vc-. From Figs. Ot. and^, we find 
that 




-LO -5 



eVL/(2j:riVp) 



5_ = 0.67 ±0.07, 



(40) 



FIG. 10: The differential conductance interference pattern 
near the lower crossing point calculated by Eq. (1321 for tun- 
neling between right movers (and similarly for left movers) us- 
ing a smooth confining potential for the upper wire, Eq. (1231 . 
Vc- — 1.4wf, Akp = 4n/L, /3 — 8. We used the numerically 
found |Af(K)|^, also shown in Fig. |H] The figure has to be 
compared to experimental Fig. |3 



similarly to the value for gi obtained by the zero-bias 
anomaly in Sec. II V Also, from Eq. (|^ it follows that 
the oscillation pattern [of the principal term G{V, B)] 
gains a tt phase shift across each suppression strip. Such 
phase shifts can also be seen in experimental Figs.|3^ and 
111. 

Finally, we compare the interference pattern predicted 
by our theory, Eq. (|32|l . with the experiment. Figs. OJi, 
0^. G{V, B) calculated using a smooth confining po- 
tential [Eq. (|23 with (3 = 8] for the upper wire is 
shown in Fig. ^| Many pronounced features observed 
experimentally-the asymmetry of the side lobes, a slow 
fall-off of the oscillation amplitude and period away from 
the principal peaks, an interference modulation along 
the V-axis, tt phase shifts at the oscillation suppression 
stripes running parallel to the field axis-are reproduced 
by the theory. 

In Fig. ^2 we repeat the calculation using (3 = 22, 
which defines potential (|23|) with a similar boundary pro- 
file near the turning points of a three-times longer wire^ 
(Here by length we mean the distance between the clas- 
sical turning points, which, as explained in Sec. IIV B II 
can be somewhat different from the lithographic length.) 
Again an agreement between the predicted fFig. lll|l and 
measured (Fig. 0} oscillation patterns is apparent. In 
Fig. ^2 a few weak side lobes also appear to the left of 
the main dispersion peaks, unlike in Fig. 1101 where they 
appear strictly to the right. In addition, the interference 
modulation in the voltage direction has sharper features 
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FIG. 11: Same as Fig.[Tnibut with Afcp = IOtt/L and P = 22, 
describing a longer junction with a similar boundary profile. 
|M(«:)|^ was correspondingly recomputed (now putting 300 
electrons per spin in the upper wire). The figure has to be 
compared to experimental Fig. |1| 



in Fig. ^2 These trends are expected for longer junctions 
as the boundaries become steeper on the scale set by the 
total length. 

Tunneling between ID channels with different Fermi 
velocities can also yield an interference modulation sim- 
ilar to that described in this section even when the 
electron-electron interactions are vanishingly small. It 
is thus important to emphasize that we suggest the spin- 
charge separation picture to explain this modulation re- 
lying on the experimental result (see Ref.|3) that the den- 
sities of modes and, therefore, the corresponding 
Fermi velocities are nearly identical. 

Using Eq. (|35|l we also studied various possible sce- 
narios when the interactions in the two wires differ. For 
example, in a situation when the upper wire is perfectly 
screened, so that T4m, = 0, there are still two veloci- 
ties present in the system, v-p and Vd, but the interference 
pattern is qualitatively very different from that shown in 
Fig. 1101 and observed experimentally [see Figs. 13 and 0]. 
Since a considerable weight of the charge-excitation con- 
tribution to the tunneling strength is shifted to velocity 
wp (which is now also the charge-excitation velocity in 
the upper wire), the oscillation pattern does not exhibit 
the pronounced vertical suppression stripes, but rather 
a much weaker modulation. The same conclusion also 
holds for intermediate regimes of relative screening in the 
two wires, when the system is not symmetric and the two 
charge-excitation velocities significantly differ. The pro- 
nounced suppression stripes are, therefore, present only if 
most of the charge-excitation tunneling weight is peaked 
at a single velocity Vc- (which is guaranteed only when 



the system is nearly symmetric). 

Taking into account 1D-2D scattering in the upper 
quantum wire will smear out the oscillation pattern by 
its convolution with a Lorentzian in the _B-direction, sim- 
ilarly to Eq. (|52|l below. The corresponding effect is, 
however, small because of the high quality of our wires, 
which have a long scattering lengthi^ 'id-2D ~ 6 /im. 



3. Upper crossing point 

In practice, since the fields necessary to reach the 
upper crossing point are quite large (e.g., 7 T for the 
<-> transition), even atomic-scale disorder in the 
junction can lead to a significant variation 5qB of the 
momentum transfer along the tunneling region. In par- 
ticular, 6qB = eB6d can be comparable with 27r/L, the 
reciprocal wave vector of the upper wire. This can sig- 
nificantly broaden the principal dispersion peaks. Fur- 
thermore, Zeeman splitting becomes about a per cent 
of the Fermi energy at these high fields and results in 
somewhat different dispersions for different spin modes. 
Away from the main peaks, however, we still expect to 
see side lobes due to stationary phases at the ends of the 
junction, similarly to the regime of low magnetic fields 
discussed above (with possibly a faster decoherence in 
the V direction than just due to the dispersion curvature 
studied in Sec. lIVB4|l . Such oscillations [with about the 
period ^] are indeed observed experimentally, as can be 
seen in Fig. |S| Because of the mentioned complications, 
we, nevertheless, do not pursue a detailed analysis of the 
conductance near the upper crossing point in this paper. 



4- Dephasing of the oscillations 

It is evident from Figs.|2t^ and0^ that the interference 
decays as is increased. A more quantitative analysis 
of this decay is shown in Figs. |3b and 2b; where the 
amplitude of the oscillations is plotted as a function of 
voltage. It is clear that the measured modulation has 
a fast-decaying envelope, which can not be explained by 
the analysis of section llVB 21 (See, for example, Eq. 
which predicts that the modulation is roughly periodic.) 

One scenario for the dephasing occurs even in the 
case of noninteracting electrons considered in Sec. lIVB II 
when we take the finite curvature of the single-particle 
dispersions into account. Let us return to the form of 
the current in Eq. Hlt)|l : 

/(X r de[\M{n+)\' + \M{n.r] . (41) 
Jo 

Correcting our previous results to take into account the 
nonlinear dispersion near the Fermi points, we now write 
K± = [fc|„ + 2me/n^]^/^ - [kli + 2m(e - eV)/h^Y/^ ± qg. 
[Using Eq. H41|) we still imply low enough bias V, so that 
the density of states in the wires are relatively constant 
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on the energy scale of elVj.] Expanding this expression 
to lowest order in curvature, we further obtain 

eV eV{eV-2e) 

[Eq. l(T^ can be recovered by neglecting the last term 
above.] The current (|41|l then becomes 



/ cx 



l'eV/2 




/ de 


M 


1 -eV/2 





eV 



eeV 



(43) 



It is easy to see now that the contribution to the 
conductance obtained by differentiating the integrand in 
Eq. (|43|) will be suppressed when the argument k of the 
tunneling matrix amplitude M{k) changes by the full pe- 
riod of oscillations An upon energy e variation between 
the integration limits ±eysup/2. We thus arrive at the 
condition for the suppression voltage Vsup- 



Ak = 



(e^sup)' 



h^v^kp 



(44) 



Approximating An « 2'k/L and translating it into the 
oscillation period in the bias direction eAV = JivfAk, 
one finally obtains 



AV 



(45) 



Using density 100 [ivoT^ for the lowest bands in the 
wires,^- we find V^up/AV^ w 7 (« 12) for the 2 /im (6 /im) 
junction. An implicit assumption in the derivation is 
that we are still close enough to the Fermi level so that 
higher-order corrections should not modify the result sig- 
nificantly [in particular, for the calculation of the matrix 
element H18|l it can still be reasonable to use the wave 
function tpuix) at the Fermi energy]. 

The result of the numerical calculation using Eq. I|43() 
and the matrix element M{k) plotted in Fig. |H1 (us- 
ing parameters characteristic for the 2 /xm sample) is 
shown in Fig. ^| Notice that when the voltage exceeds 
Vsnp ~ 7AV, so that the pattern starts dephasing due to 
the finite curvature, a beating pattern appears. It dif- 
fers from the data in several important aspects: First of 
all, the hues of suppressed G{V, B) are not equidistant. 
In addition, Kiup, corresponding to the first suppression 
stripe (on either the positive- or negative- voltage sides), 
is about twice larger than the period we observe in Fig.|3|3 
and four times larger than that in Fig.QJ), which in both 
cases is given by about 3AF. This suggests that the 
source of the beating in the experimental data is not the 
curvature of the dispersions, but rather the spin-charge 
separation mechanism discussed in Sec. IIV B 21 

Another important difference between Eq. (|41|l and the 
experiment is that the decay of the oscillations is much 
stronger in the latter. It might therefore be necessary to 
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eVL/(27:riVp) 

FIG. 12: The differential conductance interference pattern 
near the lower crossing point calculated by Eq. I|43|l . within 
the noninteracting electron picture, using the matrix element 
M(k) shown in Fig.|Hl(for /3 = 8). See text for further details. 



consider both the curvature and electron interactions in 
order to understand the fast decay of the conductance os- 
cillation amplitude with increasing voltage. Taking into 
account the curvature while bosonizing excitations of the 
interacting electronai^, leads to higher-order terms in the 
Hamiltonian. Physically this corresponds to interactions 
between bosonic excitations which therefore acquire a fi- 
nite life time. The singularities of the spectral densities 
will correspondingly be rounded, in turn smearing the 
conductance interference pattern. Further complications 
may arise from the electron backscattering which was 
entirely disregarded: While the low-energy properties of 
the system are not affected by the backscattering (apart 
from rescaling of certain parameters) since it renormal- 
izes downward in the case of repulsive interactions, the 
story at a finite energy could be different. The reason for 
this is a slow (logarithmic) renormalization flow of the 
backscattering strength. If a significant backscattering is 
present in the original Hamiltonian, it could therefore be 
still considerable at a finite energy. A detailed study of 
these effects however lies beyond this paper's scope. 



C. Zero-bias anomaly 

1. Crossing points 

It is enlightening to further study tunneling between 
ID channels at low bias when the magnetic field is tuned 
to match two Fermi-points of the wires (see Sec. IIII 
The zero-bias properties are similar near the two cross- 
ing points and, for definiteness, we choose to discuss the 
upper crossing, where the magnetic wave vector qs is 
close to k-pu + k-pi and the field changes the chirality of 
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the tunneling electrons: The tunneling is amongst the 
left movers of the upper wire and the right movers of the 
lower wire. For the ^ \li) transition, this point is 
located at i? « 7 T, see Fig. |21 The results are straight- 
forward to apply to the regime of the lower crossing point, 
as well. 

For clarity, we start by making a series of simplifying 
assumptions which will be dropped in subsequent gen- 
eralizations: First, we set the upper-wire and interwire 
interactions, V^u and Vui, to zero. Physically, this corre- 
sponds to a regime where the Coulomb interactions in the 
upper wire are perfectly screened by the 2DEG. Secondly, 
we further simplify the model by assuming a square- well 
confinement for the electronic states in the upper quan- 
tum wire and an infinitely-steep reflecting left boundary 
for the electrons in the lower wire, i.e., Uu{x) [Ui{x)] is 
constant for \x\ < L/2 [x > —L/2] and infinite other- 
wise. As we showed in the previous sections, both of the 
above assumptions are not very realistic for the purpose 
of studying the interference pattern. In the zero-bias 
anomaly regime, however, they can be a good starting 
point, at least, for pedagogical reasons. 

Electron states participating in tunneling near the 
crossing points [Eq. (^] lie close to the Fermi levels in 
both wires. It is therefore possible to calculate the cor- 
relation functions analytically using LL theory, after the 
dispersion relations in the wires are linearized. At the 
upper crossing point, we only need to retain Green func- 
tions of the left movers of the upper wire and the right 
movers of the lower wire. At zero temperature these are 
given by 



G>{x,t + iO+;x',0) 



L — >oo 



g-ifcF„Zg-r|z|/t)F 

AL sin^(z-f vpt) 

^ g-ifcF,i2g-r|z|/t)F 



27r z + vpt 
for |a;|, \x'\ < L/2, and vanishing otherwise, and 

1 



(46) 



Gf{x',0;x,t + iO+) 



27r (z — vpt)^ 



1 



{Z - Vclt) 2 



9i+S, 



z2 - {Vclt - ITcY 

91-9 



v'2 



Z 



-,12 



{Vcltf 



(47) 



for x,x' > —L/2, and vanishing otherwise, where z — 
x — x' , z' = X + x' + L, and Tc is a small distance cut- 
off. As specified above, Eq. (|46|l [Eq. (|47l) ] contains only 
the component for the left (right) movers in the upper 
(lower) wire; we have thus omitted terms proportional to 
gifeFZ^ gifcFZ ^ g-ifcFz' -vvJiich do not contribute con- 
structively to tunneling near the upper crossing point. 



The last factor in the expression for Gf is due to the 
closed boundary at x = — L/2ii2ii24i^ 

For sufficiently large voltages, eV :s> 2hvp/{giL), the 
tunneling electrons do not feel the junction boundaries on 
the time scale set by the voltage. In particular, the left 
boundary of the lower wire does not affect the dynamics 
and, effectively, electrons directly tunnel into the bulk 
of the lower wire: The last term in Eq. H47|) is close to 
unity and can, therefore, be omitted. Terms of the form 
l/^z^vtY entering Eqs. H46|l and (|47|l are dominated by 
the long-t behavior in the integral [Eq. H25I) . the voltage 
is assumed to be positive] if eV <C hin.syi{vq,T), where 
q — qB — {k-pu + k-pi). The conductance is then suppressed 
as a power law 



G{V) cx V 



(48) 



with the exponent abuik = (5/ + .9r ^ ^ This result is 
easy to generalize for the case of unscreened interactions 
in the upper wire: 



ttbulk 



E 



9,^ + 9y 



(49) 



If the interwire interactions Vui are also significant, the 
elementary excitation modes in the wires become cou- 
pled and Qfbuik has a more complicated form than that in 
Eq. (|49|) .^ Interference oscillations discussed in Sec. lIVBl 
can modulate the power-law current suppression H48|l . 
setting an upper voltage bound, eV < eAV ~ 2'Kh,v-p/L, 
for the validity of Eq. H48|l . It would therefore be 
hard to observe the exact power-law voltage dependence 
(|48|l with the exponent H49|l in the regime when eV ^ 
2hvY/{giL) (see, however, Sec. lIVC'2)l . 

If eV ^ 2hvF/{giL), electrons effectively tunnel into 
the end of the lower wire and the current suppression 
is governed by processes in the lower wire outside the 
tunneling region. In particular, details of the interactions 
in the finite upper wire do not play a role. The last term 
in Eq. H47|) now also contributes to the exponent of the 
long-i asymptotic, and a in Eq. (|48|l is given by 



O^cnd 



9r' 1 



(50) 



The upper wire, in this case, can be viewed as a point con- 
tact and the tunneling exponent is determined entirely 
by the properties of the lower wire outside the tunneling 
region. 

At a finite temperature, the time scale relevant for the 
discussion above is set by max{eV, kBT). The power law 
(|48|l should now be replaced with 



G{V, T) cx T"F„ 



eV 



(51) 



where Fa(x) is a known scaling function with properties 
Fa{Q) = const and i^^ (x) cx in the limit of a; 1 At 
low temperatures the conductance yields a low-bias dip 
extending to voltages eV ^ ksT with G{V = 0) cx T". 
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In Sec. lIVBl we showed that the conductance G{V,B) 
exhibits a characteristic interference pattern due to wave- 
function osciUations near the gates confining the tunnel- 
ing region. We can easily read out the profile of this 
pattern for the current (|25|l using the correlation func- 
tions H46|l . (|47|l in the low-energy regime considered in 
this section (namely t^z): 

where M{k) is the tunneling matrix element, Eq. (|18|l . 

The discussion in this section also holds for the lower 
crossing point, where the electrons do not change their 
chirality upon tunneling. To directly apply the above re- 
sults to this regime (for definiteness, assuming we now 
consider the transition between the right-moving elec- 
trons), we only need to redefine the distance from the 
crossing point in the field direction: q = qs + ^Fu ~ ^F/ 
(and analogously for the transition between the left 
movers). 



2. Direct tunneling from the 2DEG 

It is straightforward to generalize the main results of 
the preceding section to the regime of direct tunneling 
from the 2DEG. Eq. (|25|l stays valid in this case, but 
now is Green function for the 2DEG near the edge 
of the upper quantum well. We calculate this correlation 
function and discuss its limiting behavior at low ener- 
gies in Appendix rSl The 2DEG density of states is finite 
at the Fermi energy and, therefore, the long-i behavior 
of the one-particle Green function is G^{t) cx 1/t. If 
max(ey, /cbT) <C ?iwF^F,2D, where 7ifcF,2D is the 2DEG 
Fermi momentum and v-p is the lower of the Fermi ve- 
locities of the ID band and the 2DEG, the temperature 
and voltage dependence of the differential conductance 
are governed by the exponents H49() . with g„ = 0, or 
(|^ . depending on the relation between max(ey, /c^T) 
and 2hvp / (giL). Because in this regime we tunnel di- 
rectly from the 2DEG, interactions in the ID modes of 
the upper quantum well do not play a role, and both 
ttbuik and Qfond are determined only by the interaction 
constant gi of the lower wire. While the field dependence 
of the conductance for the direct 2DEG-lower wire tun- 
neling is different from Eq. H52(l (in particular, the con- 
ductance does not exhibit a strong oscillation pattern), 
the low-energy properties stay similar to the case of the 
ID-ID tunneling. In spite of a complicated dependence 
of G{V, B) on magnetic field, the zero-bias anomaly is 
pronounced in the data for tunneling either between dif- 
ferent ID bands or between the 2DEG and the ID bands. 

As described in Sec. IIII GI we measured the zero- 
voltage conductance dip at temperatures 0.2 < T < 2 K 
on a junction of length L — 6 /im at i? = 2.5 T. It can 
be seen in Fig. |21 that at this magnetic field, the conduc- 
tance is dominated by direct tunneling from the 2DEG, 



l^a) ^ 1^2) • Since TiVF^F, 20/^5 100 K» T, the tem- 
perature dependence of the zero-bias dip can be used to 
extract the value of the interaction constant gi for the 
band \l2)- The data points and the (best) theoretical 
fitting curves are shown in Fig. (HJ we find 

ff, = 0.59 ±0.03. (53) 

The transition point between the two lines in the plot is 
consistent with an estimate 2Tiv-p / {giLkB) ~ 0.5 K for 
the second ID mode of the lower wire, {h}- 

As a consistency check, we plot in the insets to Fig. 
curves calculated using Eq. ifCT) (taking both acnd and 
cthuik for the exponent), gi and the overall propor- 
tionality constants were independently obtained from 
the power-law temperature dependence of the bottom 
of the dip, i.e., G{V = 0,T), so that at this point we 
do not have any remaining fitting parameters. The re- 
sults show reasonable agreement with the data: When 
max(ey, fc^T) > 2hvp / {giL) the data is consistent with 
« — Oihuik while when max(ey, fc^T) < 2hv-p/{giL) it 
is more consistent with a = ctcnd- Thus, in particular, 
there is a crossover between acnd and abw/fe in the data 
for G{V) at T = 0.24 K. For voltages V - 1 meV that 
are comparable to the Fermi energies of the modes par- 
ticipating in tunneling, the power-law behavior H51|) is 
replaced by a more complex structure modulated by the 
dispersions in the wires and the upper well, see Fig. |21 

V. CONCLUSIONS 

We have presented a detailed experimental and theo- 
retical investigation of tunneling between two interacting 
quantum wires of exceptional quality fabricated at the 
cleaved edge of a GaAs/AlGaAs heterostructure. The 
study focused on revealing electron-electron interaction 
effects on the conductance interference pattern arising 
from the finite size of the tunneling region and the con- 
ductance suppression at low voltage. 

In the analysis of the data the finiteness of the junc- 
tion plays a central role. Breaking translational invari- 
ance, the boundaries give rise to secondary dispersion 
peaks in the dependence of the conductance on voltage 
bias and magnetic field. Smooth gate potentials result 
in a strongly asymmetric interference profile, while the 
Coulomb repulsion in the wires leads to spin-charge sep- 
aration which, in turn, modulates the conductance oscil- 
lation amplitude as a function of voltage bias. 

An interplay between the electron correlations in the 
wires and the finiteness of the junction length also results 
in different regimes of the zero-bias anomaly. At the low- 
est voltages, the upper wire is effectively a point-contact 
source for injecting electrons into the semi-infinite lower 
wire. On the other hand, at higher voltages, electrons ef- 
fectively tunnel between the bulks of the two wires along 
the length of the junction. 

Using the temperature dependence of the zero-bias dip, 
we found the value of the interaction parameter gi = 
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vpi/vci for band IZ2) in the lower wire to be 0.59 ± 0.03. 
From the ratio between the long (due to spin-charge 
separation) and slow (due to upper-wire confinement) 
scales of the conductance oscillations, we also extracted 
the interaction parameter g- — v-p jvc- corresponding to 
the antisymmetric charge-excitation mode in the lowest 
bands and |^i) of the biwire to be 0.67 ± 0.07. 

While (7_ and gi have similar numerical values, these 
quantities should be contrasted: gi is the interaction pa- 
rameter of the channel IZ2) in the lower wire, which 
is screened by other ID states in the wires as well as the 
2DEG of the upper quantum well. on the other hand, 
is a parameter characterizing the (antisymmetric) charge 
mode in the coupled and |^i) channels of the two 
wires, which is relatively weakly screened by the 2DEG 
since the latter has a smaller Fermi velocity (being, never- 
theless, still larger than the Fermi velocity of This 
can explain why g_ and gi are comparable while [^2) has 
about half the Fermi velocity of and [The in- 
terwire interaction would only enhance the mismatch as 
it reduces Wc-, see Eq. 

Similar values for the interaction parameter in the 
range between 0.66 and 0.82, were found in Ref. |23 for 
single cleaved-edge quantum wires by measuring the tem- 
perature dependence of the line width of resonant tun- 
neling through a localized impurity state. Spectral prop- 
erties of the same double- wire structure as reported here 
were investigated in Ref. |^ also indicating comparable 
values of 5, about 0.75, for various intermode transi- 
tions. An interaction parameter g k, 0.4 was found for 
GaAs quantum-wire stacks in resonant Raman scatter- 
ing experiments;-^ the smaller value of g there can be 
attributed to much lower electron densities and no screen- 
ing by the 2DEG as in our measurements. 
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and A. Stern. This work was supported in part by the 
US-Israel BSF, NSF Grant DMR 02-33773, and by a re- 
search grant from the Fusfeld Research Fund. OMA is 
supported by a grant from the Israeli Ministry of Science. 



APPENDIX A: INDEPENDENT-MODE 
APPROXIMATION 

In our analysis we treat different ID bands in the 
wires as independent and disregard interband interac- 
tions. While this is a convenient approximation for theo- 
retical investigations that has been often assumed in pre- 
vious workS)SiS*^ it needs to be further justified. Tunnel- 
ing into multimode ID wires was considered in Ref. |23, 
It was shown that low-energy tunneling into the edge of 
a semi-infinite wire with A'^ bands is governed by the tun- 
neling density of states exponents ai such that the differ- 
ential conductance (at zero temperature and low voltage 
y) is given by G oc X]i=i where ti is the tunnel- 

ing amplitude for the zth mode. In the independent-mode 
approximation with interactions described by Hamilto- 
nian for each mode, these exponents are given by 



Eq. H50|l with the parameter g describing interactions in 
each mode. On the other hand, in a more realistic picture 
one deals with an interaction Hamiltonian 



N 

E 



dxpi{x)pj{x) (Al) 



which takes into account the interband coupling. Here, 
Vq is the zero-momentum Fourier component of the in- 
teraction potential V{x) — Vi:)5{x) and pi is the electron 
density in the ith band. The exact form of the potential 
is not important as we are only interested in the long- 
wave-length quantum fiuctuationsi2^ 

The exponents are given bjiSS Ui = (X]/!=i la^i/'^i) ~ 
where Vi is the Fermi velocity of the noninteracting ID 
electron gas at the density of the ith mode, s; is the veloc- 
ity of the Zth soundlike excitation in the presence of the 
potential V{x), and 7^; characterizes coupling between 
the ith and /th noninteracting modes after the interac- 
tion potential V{x) is switched on. In the case of a single 
transverse mode with spin degeneracy, A^ = 2, 7^^^ = 1/2, 
and si — wp v/l + 2 Vo / {nhvp ) , S2 = wf are the charge- 
and spin-excitation velocities, respectively. For a general 
N, the velocities si are given by roots of the equation 



«2 _ ,,2 = 



-1 'I 



(A2) 



and the coefficients are given by 



ll 



2\2 



N 

E 



(A3) 



In our systemj^ the Fermi velocities of the highest oc- 
cupied bands are very different (e.g., the highest trans- 
verse mode has twice the velocity of the next lower- 
lying mode). Furthermore, since the interaction Vq ^ 
max(?iw.i) is not too large, the correction to the expo- 
nents ai due to the interband coupling is expected to be 
relatively small. One can accommodate for this correc- 
tion by slightly renormalizing the interaction constants 
g, viewing it as a mutual interband screeningi^^ 

Also, it is safe to disregard intermode transitions as 
they are determined by the Fourier components of the 
interaction with a large wave vector fc ~ fcp, which 
are small for a smooth long-range potentiali^ The weak 
backscattering within each spin-degenerate mode can be 
further renormalized downward at low energies in the 
physical case of repulsive interactionsii^ 



APPENDIX B: DIRECT TUNNELING FROM 
THE 2DEG 

In order to describe the V and B dependence of the 
conductance for direct 2DEG-lower wire tunneling, we 
approximate the Green function of the top quantum well 
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by the edge Green function of a 2D electron gas occupying 
a half plane y > with x extended from — cxd to oo. 
We assume the potential is V{x,y) — for y > and 
V{x, y) — oo for y < 0. Therefore, we find 

iG>(x,y,t;x',y',Q) ^ / dpe^P^^"^'' / dk 

I" J-oo Jo 
X sin(fcy) sin(fcj/')0(e)e~"*^'* , 

(Bl) 

where e = h^{p^ + fc^ — kp)/{2m) is the energy and fc| is 
the Fermi wave vector of the 2DEG. 8(e) is the Heaviside 
step function. When we calculate the tunneling current, 
y and y' run from to ^, the width of the tunnel junc- 
tion (i.e., the extent of the ID mode of the lower wire in 
the direction perpendicular to the cleaved edge). We set 
{y,y') £,/2 and approximate sin(fc^/2) w k£^/2 assum- 
ing fep < 1/^- In the frequency domain, Green function 
G>{z,uj) = J^^dte''^*G>{z,t), with z ^ x - x' , then 



becomes 

iG>(z,w) = ^/ dpe'P^ dkk^S{e-uj)eiuj) . 
27r?i J_oo Jo 

(B2) 

In the limit of small positive frequencies it reduces to 

iG>(z,w^O+) = m-^ / ^ dpe'P'Jkl-p^ 
2ttTi J^kj, V 

2 Ji{kFz) 

= m{£,kF) -7rr-r , (B3) 

InkYZ 

where Ji is the first-order Bessel function of the first kind. 
In particular, since Ji{x) oc x when x — *■ 0, the density 
of states is finite at the Fermi energy and G^{t) cx 1/t 
as t oo. Furthermore, from the low-energy form of 
the 2DEG Green function [Eq. |jB2Il] it follows that the 
relevant range of z in integral is l/fcp rather than 
1/ max(g, F/wp) as in the case of the ID-ID tunneling. 
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